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Section  I 
Introduction 

For  the  linearized  acceleration  potential,  as  well  as  the  velocity  potential  equation  in 
three-dimensional  unsteady  flows,  there  exists  a  fundamental  solution  in  a  closed  form. 
This  makes  it  possible  to  formulate  the  problem  of  the  flow  over  a  wing  (or  also  tor  the 
whole  airplane)  in  terms  of  an  integral  equation,  by  which  the  problem  is  reduced  from  three 
space  dimensions  and  time  to  two  space  dimensions  and  time.  Under  special  conditions,  for 
instance  the  flow  over  an  oscillating  wing,  the  problem  has  been  solved  before;  an  example 
is  the  doublet-lattice  method.  In  these  approaches,  valid  for  steady  or  oscillatory  flows,  the 
integral  equation  formulation  does  not  appear  explicitly.  It  is  by-passed  by  an  immediate 
discretization  of  the  problem  combined  with  a  specialization  of  the  fundamental  solution. 
By  this  procedure  one  avoids  the  need  to  deal  with  the  singularities  which  by  necessity 
occur  in  the  integral  equation. 

The  authors  of  Reference  1  (Guderley  and  Blair)  believe  that  for  general  applications 
which  deal  with  the  problem  in  the  time  domain,  the  explicit  integral  equation  may  be 
advantageous.  However,  to  bring  the  integral  equation  into  a  form  suitable  for  a  numerical 
approach  it  is  necessary  to  carry  out  transformations  which  weaken  certain  singularities 
otherwise  encountered.  The  analytical  work  carried  out  for  this  purpose  in  Reference  1  is 
rather  formidable.  By  performing  the  necessary  steps  in  a  different  order,  Professor  Marc 
H.  Williams  has  accomplished  considerable  simplifications.  Moreover,  Professor  Williams' 
approach  leads  to  a  simplification  of  one  of  the  final  equations  and  makes  a  complete 
exploration  of  the  relations  between  the  steady  and  the  unsteady  problem  possible.  In 
addition,  it  opens  the  way  to  a  rather  intriguing  development,  namely  the  treatment  of 
the  integral  equation  by  means  of  Lorentz  coordinates. 

The  present  report,  prepared  by  K.G.  Guderley  and  Mark  R.  Lee,  follows  the  notes 
kindly  sent  to  Guderley  by  Professor  Williams.  Some  of  the  intermediate  steps  have  been 
filled  in.  to  facilitate  the  study  of  the  report  for  a  reader,  who  on  one  hand  does  not  want 
to  accept  the  results  on  faith  but  on  the  other  hand  has  only  limited  experience  in,  or  time 
for.  the  necessary  mathematical  manipulations. 


I 


Section  II 

Basic  Equations  in  Acoustic  Coordinates 

Let  x  .  y  .  ~  bo  a  Cartesian  system  of  coordinates  which  is  at  rest  with  respect 
to  the  undisturbed  air.  The  superscript  zero  is  applied  because  at  some  later  stage  the 
coordinates  x.  y.  z  will  be  introduced  for  a  Cartesian  system  fixed  with  respect  to  the  wing. 
With  respect  to  this  (no  superscript)  system,  the  undisturbed  air  moves  from  left  to  right 
with  constant  speed  in  the  x-direction.  The  planform  of  the  wing  lies  in  the  x,y-plane.  Let 
i.j.k  be  unit  vectors  in  the  x  .  y  .  z'  or  x.  y.  z-directions  respectively.  Let 


x  =  x  i  +  y~ ]  +  z°k 


(2.1) 


Accordingly  x  represents  a  vector  with  components  xc ,  y° ,  zc .  The  reader  should  not 
confuse  the  vector  x  with  its  component  in  the  i  direction,  denoted  by  x  .  In  the  x  .  y" .  z"- 
system  the  wing  moves  with  constant  velocity  in  the  direction  of  decreasing  xc .  In  this 
process  it  sweeps  out  a  strip  in  the  x~,yc  plane  which  has  a  width  (i.e.  an  extension 
in  the  y  -direction)  equal  to  the  span.  In  the  xc-direction  this  strip  extends  from  the 
position  of  the  wing  at  the  time  when  the  unsteady  motion  started  to  its  position  at  the 
current  time.  In  the  present  analysis  viscosity  is  disregarded  and  the  air  is  considered  as 
isentropic.  Since  the  initial  state  of  the  air  is  free  of  vorticity,  the  Helmholtz  vorticity  laws 
are  applicable,  and  the  velocities  can  be  described  as  the  gradient  of  a  potential  (denoted 
by  <2>).  The  equation  for  the  potential  is  obtained  by  the  requirement  that  the  equation  of 
conservation  of  mass  be  satisfied.  We  make  the  further  assumption  that  the  perturbation 
generated  in  the  air  (which  was  originally  at  rest)  is  small.  Then  the  equation  for  the 
potential  reads 


a2,  V  2  6 


d2<t>z 


(2.2) 


Here  a  is  the  velocity  of  sound  in  the  unperturbed  state,  and  V°2  is  the  Laplace  operator 
in  the  x  .  y  .  z  -system.  The  del  operator  ?°  is  defined  by 


V° 


.  d  .  d  ,  d 
ldx  +Jdy  +  kbz< 


Let  f  lie  a  fixed  vector  in  the  x  .  y  .  2  -space. 

£'  =  +  jr]  +  k( 


a  nd 

R  =  R  =  [(x  -  n2  +  (y  -  n  )2  +  (r  -  cr)2]’  2 


(2.3) 


(2.4) 

(2.6) 


2 


(2.7) 


where  the  prime  denotes  the  differentiation  with  respect  to  the  argument.  For  R0  — *  0, 
one  has  r  —  t.  and  the  last  expression  reduces  to  h[t).  Accordingly,  h(t)  is  the  mass  flow 
per  second.  This  is  the  interpretation  of  the  function  h(t).  While  this  discussion  has  no 
bearing  on  the  following  development,  it  does  explain  the  factor 

At  the  point  given  by  x'  —  the  potential  equation  is  not  satisfied,  for  there  is  a 
mass  flow  emerging  from  this  point.  The  flow  will  be  represented  by  a  superposition  of 
particular  solutions.  Eq.(2.7),  with  centers  at  different  points  The  functions  h°  are 
allowed  to  be  different  for  each  value  of  £c.  This  is  expressed  by  w’riting  h0(tc,r).  Notice 
that  £  also  occurs  in  the  definition  of  r,  thus  appears  in  h  twice,  implicitly  (in  r)  and 
explicitly  as  first  argument.  The  singular  points  x  =  £G  cannot  lie  inside  of  the  flow  field; 
and  they  can  appear  only  wdthin  the  strip  of  the  x' ,y°-plane  swept  out  by  the  wing.  The 
component  of  the  vector  is  therefore  zero.  This  assumes  the  wake  to  be  flat. 

The  unsteady  up  and  dowm  motion  of  the  wing  generates  at  its  surface  a  velocity  in 
the  z  -direction  which  is  the  same  on  the  upper  and  lower  side  of  the  wing.  The  field  of 
the  s'  components  of  the  velocity  is  therefore  symmetric  with  respect  to  the  x°,y°-plane. 
The  potential  is  then  anti-symmetric  with  respect  to  this  plane.  This  implies  for  points 
in  the  x  .  y  -plane  outside  of  the  strip  swept  out  by  the  wing  that  the  potential  is  zero. 
The  source  potential  of  Eq.  (2.7)  is  spherically  symmetric:  therefore,  it  is  symmetric  with 
respect  to  the  x  .y  -plane  if  £  =  +  jy  .  An  antisymmetric  expression  is  obtained 

by  a  differentiation  with  respect  to  s'  .  The  resulting  expression  can  be  interpreted  as  a 
doublet  (of  moment  -1)  oriented  in  the  positive  .?  direction  (actually  this  interpretation 


is  unimportant).  Thus  we  represent  the  potential  by  a  superposition  of  “doublets’" 


6d(x-.£-J) 


-1  8  h  (s  •’) 
4  7T  dz  R 


where  £'  =  ;£  +  jr)°,  is  restricted  to  the  above  mentioned  strip  of  the  x°,yc-plane  for  the 
remainder  of  the  report.  Accordingly,  the  potential  will  appear  in  the  form 

=  -iiSF  J  <*•»> 

The  uinbral  variables  are  £"  and  y' .  One  must  remember  that  £  occurs  implicitly  in  r. 
The  integration  must  be  extended  over  the  above  mentioned  strip.  So  far,  the  function 
/r(£:,r)  is  unknown.  One  has,  because  =  0 

dRo  _ 

8z°  R0 


therefore, 


» i.-.//  -- 


hcU°,r) 


R3o 


hr{(0.T) 

a0Ro 


d(cdyc 


(2.10) 


where  h °  denotes  the  derivative  of  hz  with  respect  to  the  second  argument  r.  In  potential 
theory  one  determines,  for  a  sheet  of  sources,  the  normal  component  of  the  potential  in 
terms  of  the  local  source  intensity.  The  expression  corresponding  for  the  sources  in  the 
present  case  is  given  by  Equation  (2.7).  Equation  (2.9)  is  then  the  expression  for  the 
component  of  the  potential  gradient  in  the  z-direction.  We  shall  show  by  applying  an 
analogous  procedure  that 


lim  d>°(t'x +  jy  +  kz~ ,t)  -  lim  <t>°(ixr  +  jy°  +  kz° ,t)  =  hz(ix°  +  jy°,t)  (2-11) 
2'---0  2°— -o 

Here  we  repeat  the  main  steps  of  the  argument.  Consider  a  point  given  by  xc  =  ix°  + 
jy'  +  kzn  and  determine  the  limit  of  Eq.(2.10)  as  z°  — *  +0  while  x°  and  yc  are  fixed.  We 
divide  the  region  of  integration  into  an  inner  and  outer  region  by  a  small  circle  with  radius 
p  around  the  point  £  =  x  .  y  =  y  .  Outside  the  circle  R0  ^  0;  the  integral  (over  the 
outside  region)  vanishes  in  the  limit  z(  — > ►  0.  To  carry  out  the  integration  within  the  small 
circle  we  introduce  polar  coordinates. 


£  —  x  —  p  cos  0 
y  -  y  —  p  sin  0 


4 


The  dominant  term  of  the  integrand  for  small  p(  is  h 
The  contribution  of  this  term  is 


where  £  =  ix  +  jy  ■ 


47 


The  integral  is  rewritten.  Setting  p/r  -  p.  one  obtains 


r  ~  p  ,  ['p  P: 

Jo  ( p 2  +  ;  2)3  2  P  ~  Jo 


( P 2  +  1 


3  2 


dp  =  -(p2  +  1) 


nr1  4’ 


Hence 


lira  ( 
-OJo 


P  Z  p 

0  (p2  +  ^2)^d(>  =  1 


VN  ith  Equation  (2.9).  and  the  discussion  following  it.  one  now  obtains 


2:  --0 

Since  the  potential  is  antisymmetric  one  has 


lim  o  (ix'  +  jy  +  kz°.t)  =  0hc{ix'  +  jy\t) 


lim  d  (ix'  +  jy'  +  kz\t)  =  -  h°(ix:  +  jy°,t) 
’  ‘  ■*  0  £ 


Accordingly  h  (ix  +  jy  .  <)  represents  the  difference  of  the  potential  between  the  upper 
and  the  lower  side  of  the  doublet  street.  at  the  point  given  by  x°  =  ( ix°  +  jy°).  The  same 
result  can  be  obtained  by  a  simpler  argument.  According  to  the  discussion  leading  to 
Eq.(2.9)  the  function  o  (x  .t)  can  be  interpreted  as  the  2C  component  of  the  potential 
generated  by  a  source  sheet  of  strength  hc(ix°  4-  jy  ,/.).  We  showed,  following  Eq.(2.8), 
that  the  mass  flow  per  unit  time  of  a  single  source  at  a  point  is  h°(£\t).  The  mass 
flow  per  unit  time  emanating  from  an  element  dA  surrounding  the  point  =  ix°  4-  jyc 
is  then  h°(d  ,t)dA.  The  mass  flow  is  also  expressed  by  the  difference  of  the  zc  velocity  at 
the  upper  and  lower  sides,  multiplied  by  dA.  But  the  velocity  components  are  given  by 
lirn^d  (ix  +  jy  +  kz  .1)  hence  the  result  given  in  Eq.(2.11). 

The  motion  of  the  wing  has  a  2  component  of  the  velocity,  which  will  be  denoted 
w ( x  .y  .t).  In  the  linearized  approximation  this  is  the  velocity  component  normal  to  the 
wing  and  as  such  it  is  transmitted  to  the  air.  Therefore 


rr(x  .  y  .1)  =  lim 


ft 


dd  (x  ,  /) 
dz 


(2.12) 


where  «c  is  known  on  the  wing  surface 


do  [j  .  I) 
dz 


1  d2 
47 :  dz  2  J 


1 1,  /?  'r'CWfW 


(213) 


T  ~  t- 


R 


R 


x  — 


~  [(-c  -  $  )2  +  (y  -  »/  )2  +  c2 


Furti 


H'rinore 


dR 

dz 

()T 

dz 


R 


a  R 


Carrying  out  the  differentiation  with  respect  to  z  we  obtain 

--Iff 


do  (  /  .  / ) 
dz 


d2  1 
dz  2  R 


h"(Z''r)+  a  ra  Mf  ,r)  +  a2R?h}T^: 'T)) 


di'dr]c 

2.14) 


where 


£  -  »T  +  jn 

According  to  Equation  (2.12)  one  must  carry  out  the  limiting  process  z’’  — *  0.  If  the  point 
( x  .  y  )  lies  outside  the  region  of  integration,  this  can  be  done  immediately  and  one  obtains 


do  {Jr  -t) 
dz 


z  o 


R 


\h  (£.,)  +  -A-hUf'-T) 

■  a,  Rr 


d£r  dr] 


(2.15) 


This  form  of  the  equation  will  be  useful  in  conjunction  with  the  acceleration  potential. 
There  one  encounters  the  task  of  evaluating  this  expression  for  points  of  the  x.y-plane 
upstream  of  the  wing. 

For  £  =  ix  +  jy  .  one  obtains  R  .  =  0  and  then  it  is  not  possible  to  carry  out  the 

limiting  process  c  —  0  lirectly.  Actually  some  contributions  to  the  integral  cancel  each 
other.  This  is  shown  by  .  ns  of  certain  transformations  (in  essence  they  are  integrations 
by  parts).  We  present  here  the  approach  of  Professor  Williams.  He  points  out  that 
satisfies  the  Laplace  equation  in  the  x  .  y  .z  system,  therefore; 


dl  /  1  \  /  d2  /  1  \ 

dz  2  \  R  I  ~  \  ds  2  \  R  I  + 


d2 

dy  2 


1 

R 


but  as 


R 


+  (y  -  n  )2  +  (.:  )2  2 
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One  also  has 


with 

—  i,  .  B  .8 
-  l8i-  +3  drr 

The  superscript  o  indicates  that  the  del  operator  is  defined  in  the  f  ~ ,  y  -system  (as  opposed 
to  £.  p-systern);  the  subscript  2  emphasizes  that  it  is  defined  in  the  two-dimensional  £c: .  p  ,- 
space  (rather  than  in  the  £ \  rj  .  ^  -space).  The  first  term  in  the  integrand  of  Equation 
(2.14)  therefore  assumes  the  form 


h'(C- 


B2 

]  8  zT‘ 


This  is  rewritten  by  means  of  the  following  relation,  easily  verified  if  one  writes  it  out  in 
components.  Let  a(£°,7?c)  be  a  scalar  field  and  u(£c,t?°)  a  vector  field.  Then  one  has 


aVC2  ■  v  =  Vj  •  (a?)  —  •  v 


Accordingly 


h'ir.r) 


In  forming  ^2^' (£'■”)  one  remember  that  the  vector  £c  occurs  in  h°  as  the  first 

R 

argument  and  also  implicitly  in  r  =  t - To  make  clear  which  partial  derivatives  are 

meant,  one  best  writes  h'  (^  .r)  =  h  (£°,7?',r)  and  applies  the  usual  notation  for  partial 
derivatives.  Then 

V2(Msf  ,  t))  =  ih\  +  jhr,  +  hTVr2r 


(even  for  zc  ^  0).  L°t 

'T(^wr)  =  ih\,  +  jh'  o 

and 

P  -  i(C  ~  )  +  j{nr-  -  y~) 


then  even  for  z  /  0 


V,R 


P 


R 


(2.16) 


(2.17) 
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±)- 


(2.18) 


ami  by  E<iuation  (2.8) 


V2r  =  - 


Therefore.  with  />  ■  p  —  \p  \~  =  [R*  —  z  ’*) 


ih  (  1 


V2  •  [h  (£  -r)V2  ))  -3  •  -?U  •")  +  hrU^T)  + 


-tituting  thb  expression  into  Equation  (2.14)  one  obtains: 


+  hr^:-T)  +  alR^7^ 'T)\  dedr}° 

Notice  that  these  transformations  hold  for  zc  ±  0.  Thus  R  /  0  and  the  expression  on  the 
left  in  the  following  equation  can  be  transformed  by  the  Gaussian  integral  theorem. 


//,  -<v[*  (*)]«’*•  =  -//‘(eaFi (£)  •* 


(2.19) 


where  r>  is  the  unit  vector  in  the  direction  of  the  outer  normal  to  the  boundary  of  the  area 
of  integration,  ds  is  the  line  element  and  5  is  the  contour  of  the  region  A0.  But  we  found 
that  outside  of  the  area  of  integration  o  (£° A)  =  0  if  lies  in  the  xL\y°  plane  and  that 
h  is  the  jump  of  o  at  the  plane  e  =  0:  there  is  no  jump  of  0°  outside  the  strip  of  the 
i  .y  plane  swept  out  by  the  wing.  Moreover,  0°  is  continuous,  as  one  moves  across  the 
boundary  of  the  region  .4  (otherwise  there  would  arise  a  delta  function  in  the  velocity). 
Therefore  h  ( f  .  /  )  =  0  at  the  boundary  of  Ac  and  the  expression  shown  in  Equation  (2.19) 
vanishes.  This  procedure  eliminates  the  worst  of  the  singularities.  Equation  (2.12)  is  now 


written  in  the  form 


where 


■(7  ./)  =  wt,(xrj)  +  wJx  ./) 


=  4  U  IA 


(2.20) 


(2.21) 


h-  [f  .  I)  -  I 


L  /  jA  +  «- -r)]rf€rf,  (2.22) 
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provided  that  h  (£' .  r)  e  0  at  the  boundary  of  the  region  .4,.  The  expression  q  =  ih}c 
j’/t  is  cefined  in  Equation  (2.16)  and  can  be  interpreted  as  follows.  We  found  that 
o  ( £  . //  .  r)  =  (£  .  r/  .  r).  therefore  ^  c ;  +  )  gives  the  velocity  component  on  the 

surface  at  the  upper  side  of  the  sheet  and  —  £( //?  c  -f  jh '■)  is  the  velocity  at  the  lower 
side  The  strip  on  which  h  j-  0  is  therefore  a  shear  layer,  which  can  be  interpreted  as  a 
vortex  sheet  (Hence  the  subscript  v).  The  two  velocity  components  have  opposite  sign.  The 
vorticity  vector  lies  in  the  £  .  r)  plane  and  is  perpendicular  to  the  vector  +  jh'^ ).  The 
p  •  (  f '  .  T ) 

expression  ' ^  —  can  be  obtained  by  the  Biot-Savart  law.  Not  much  is  gained  by  this 
interpretation.  An  isolated  vortex  element  does  not  exist,  because  it  violates  div  curl  v  — 
0.  An  infinitesimal  vortex  ring,  according  to  Stokes'  vorticitv  law,  is  equivalent  to  a  doublet 
distribution  over  a  surface  which  has  the  vortex  ring  as  boundary,  and  this  brings  us  back 
to  the  starting  point  of  these  investigations. 

For  2'  —  0,  and  £  —  x° .  the  distance  Rc  — >  0  and  one  still  encounters  a  singularity 
because  of  the  denominator  R?  in  Equation  (2.21).  We  rewrite 

«v(7.<)  =  i‘;//5UV).j|4W  = 

+jy,‘)f  f  'gid£  di' +  fi  J  /(t(C-r)  - 5(<v  +  jy'-T))- (223l 


In  the  second  term  the  singularity  vanishes.  This  becomes  evident  if  the  integration  is 
carried  out  in  polar  coordinates: 

—  x  =  pros  6 

r)  —  y  =  psin  6 
d£  dr)  —  pdpdO 


where 


P  =  \p\ 


Now  ("T(£  .  r)  -  ^[ix  +  jy  ./))  is  O(p).  We  therefore  obtain  an  integrand 


>2  +  *2)3’2 


which,  in  the  limit  z  — •  0.  becomes  O{l)d0.  In  the  first  integral  in  Equation  (2.23),  the 

introduction  of  polar  coordinates  would  give  ^  ,  dpdO.  Carrying  out  the  limiting 

(P2  +  ~  2)  ’  1 

process  r  —  0  at  this  stage,  one  obtains  O  |  j  dp.  which  does  not  converge  for  p  =  0. 
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For  ;  7^  0  this  integral  is  rewritten  using  Equation  (2.18)  and  subsequently  the  gradient 
theorem. 

fj  tv  =  ~  j =  f  knd* 

where  rl  is  the  unit  vector  m  the  direction  of  the  outer  normal  and  ds  is  the  arc  length 
of  the  perimeter  of  the  region  under  consideration.  At  the  boundary  curve,  along  which 
the  integration  on  the  right  is  carried  out,  i?,  7^  0.  so  that  one  can  perform  the  limiting 
process  2  —  0.  then  R  —  |p[.  The  gradient  theorem  is  readily  derived  by  carrying  out 
tlie  integrations  in  the  following  expression 

/  /  i*  dx  +  dxdy  —  j  u(idy  —  jdx)  —  j  unds 

In  evaluating  ?r;j.  Equation  (2.22),  in  the  limit  zc  — >  0,  a  technique  familiar  from  potential 
theory  is  app!  d  again.  Consider  a  point  given  by 

xc  =  x  i  +  y'  j  +  zck 


By  drawing  a  small  circle  with  radius  p  around  the  point  =  x° .  =  yz  of  the 

piane.  we  divide  the  region  of  integration  into  an  inner  and  an  outer  region.  In  the  outer 
region  R  0.  and  the  limiting  process  z L  — >  0  gives  zero  immediately.  If  pc  is  sufficiently 
small  and  for  2"  — »  0.  h:T (£'- ,  r)  can  be  replaced  by  h? (ix°  +  jy°,  t)  =  h^ix0  4-  jy°  , /).  Thus 
one  obtains  as  contribution  of  the  first  term  in  Equation  (2.22)  for  1 up 

:  2 


~^ix  j  a7Ridiiri 


In  polar  coordinates  w-ith  polar  radius  p,  one  obtains 


ht(x.yj)  [»°  z'2t 


_  hf  (x.y.t)  C 

ak  Jo 


( p 2  +  Z^2 


d-P 


Hence  with  [pjz  )2  =  u 

ht(x.y.t)  fp 


r 

Jo 


z2p 


(pi  _|_  2&2)2 

ht{x.y.t) 


dp  =  - 


ht{x,y.()  /•“ 


2a 


r 

Ju^O 


*'2  dv 

(u  +  lji 


«  i :  i2 


2»  ((»"+l))).  0 


iim 


h,(x.y.t 


2ira 


Ji 


t>  z  1  pdOdp  _  ht(ix  +  jy  .t) 
(p2  +  z'i)  2  2a 


Then 


Treating  the  second  term  in  Equation  (2.22)  in  a  similar  manner:  one  obtains  zero.  Thus 


u’p(x  J)  =  - 


2<i 


For  an  interpretation  of  this  term  wo  observe  the  following.  If  w(x' .  y: .  +0. /)  does  not 

depend  upon  .r  and  u  thru  the  same  holds  for  h(x  .  y  . /).  The  integral  equation  then 
1  (lh  h 

reduces  to  »■(/)  =  9-  ^  .  As  ?/-(/)  =  o.  (l).  and  ^  =  o  (x  .y  .+()./).  one  then  obtains 

o.  —  —  -  o(  .  This  relation  holds  for  a  plane  wave  propagating  in  the  z:'  direction.  Tins 
is  the  approximation  of  piston  theory  (therefore  the  subscript  p). 

In  the  acoustical  approximation,  the  perturbation  pressures  are  expressed  (from  the 
Bernoulli  equation  for  unsteady  flow)  by 


A p  =  ~PocOt  =  h{(Zc 

where  px  denotes  the  density  of  the  undisturbed  air.  The  pressure  difference  between  the 
upper  and  the  lower  side  of  the  wake  is  0.  Therefore  ht  =  0  at  points  of  the  wake.  In  the 
x  .  y  .z '  system  the  value  of  hc  remains  constant  for  points  of  the  wake.  It  is  determined 
by  the  value  it  assumed  at  the  moment  when  the  trailing  edge  swept  over  this  point. 

In  the  present  formulation  it  has  been  assumed  that  the  wing  moves  within  the  xz ,  yz- 
plane.  It  need  not  move  in  a  straight  line  or  with  constant  velocity.  For  such  a  situation, 
or  even  if  the  wing  motion  is  not  restricted  to  this  plane,  the  x  ,y  .z;- system  is  the  ap¬ 
propriate  system  of  coordinates.  (It  is.  of  course  assumed  that  the  motions  are  compatible 
with  linearized  theory.)  In  the  more  general  cases  one  will  have  a  wake,  which  does  not 
lie  in  the  x  .y  plane  but  the  condition  —  0  for  points  of  the  wake  still  applies.  If 
the  wing  does  not  lie  within  the  x.y  plane,  then  one  must  introduce  doublets  oriented 
in  the  direction  normal  to  the  wingfrather  than  in  the  r  direction).  The  wing  motion 
and  deformation  always  determines  the  velocity  component  in  the  direction  normal  to  the 


1 ) 


wing. 


Section  III 
Wing  Coordinates 

For  a  wing  moving  with  constant  speed  J'  m  tin'  x-direction.  it  is  convenient  to  trans¬ 
form  the  equations  to  a  Cartesian  system  fixed  with  respect  to  the  wing.  Accordingly  we 

s  =  .r  I' I 

y  =  y 

/  =  / 

The  manner  in  which  £  (and  //  )  is  transformed  into  wing  coordinates  is  based  on  the 
following  reflection.  The  function  h  appears  in  the  integral  equation  with  a  third  argument 
r.  We  want  a  newly  defined  function  h(£.ri.r)  to  refer  to  a  fixed  location  on  the  wing. 
Then  the  relation  between  £  and  £  is  determined  by  the  location  of  the  wing  with  respect 
to  the  £  .  i]  system  at  the  time  given  by  r.  Therefore  we  have 

C  —  C  ,  T’~ 

S  —  S  T  <  I 

V  =  V 

M£-  >1-  r)  =  h  (  £  —  It.  //.  r  j 

As  far  a'  the  integral  is  concerned  this  is  a  transformation  from  the  variables  £  and  ry 
to  £  and  iy  In  this  transformation  x  .  y  .:  and  /  play  the  role  of  parameters,  they  enter 
through  the  definition  of  the  auxiliary  function  r. 

T  =  t  -  R  /«  (3.1) 

with 

R  -  [(£  -t  \2  +  {ti  -  y  )l  4  -r  2]‘  (3.2) 

In  at  least  one  expression  it  js  simpler  to  form  the  limit  ;  —  0  after  the  transformation 

ha~  been  carri<  d  <ei:  Tin  :  rnti'format  inn  from  £  ./  r<m  .//  i' straight  forward .  The  inverse 
transformation  i-  le"  so.  W  <  introduce 


M  -  I  i,i 


l" 


then 

{C-X  °)  =  {l-x)+U(t-T) 

Hence  with  Equation  (3.1) 

(f  ~*c)  =  U-x)  +  MR,_ 

substituting  Equation  (3.3)  into  Equation  (3.2)  one  obtains 

Rl  =  ((£  -  x)  +  MR 0)2  +  (rj  -  y )2  + 

This  is  a  quadratic  equation  for  R0.  Writing  it  in  detail  one  has 

Rl(l  -  A/2)  -  2 M{£  -  x)R0  =  (f  -  x)2  +  (r?  -  y)2  +  z2 
The  result  is  written  in  the  following  form.  Let 

Q2  =  1  -  M2 


(3.3) 


(3.4) 


and 

Rc  =  [(£  -  x)2  +  $2(rj  -  y)2  +  P2z 2] 2  (3.5) 

then 

Ro  =  ~  [M(e  -  X)  +  Rc]  (3.6) 

From  now  on  R0  and  also  r  are  considered  as  functions  of  x.y,f,77,  and  t.  Substitut¬ 
ing  Equation  (3.6)  into  equation  (3.1),  one  obtains  r;  subsequentally  one  finds  f°.  The 
transformation  of  r?°  and  z°  are  trivial.  According  to  Equation  (2.17) 

P  =  *(£°  -  x°)  +  j(rj°  -  y°) 

Then  with  Equation  (3.3)  in  terms  of  the  new  coordinates 

p  =  t(f  -  x  +  MR0)  +j{ri  -  y)  (3.7) 


Substituting  one  obtains 


P  =  02  ( ^  ~  x  +  MRe)  +  ■?(»?  -  y) 


One  obtains  from  Equations  (3.5)  and  (3.6) 


dR0  _  1  1  ,  . 

di  ~  p  rc  \MRr  +  * 


(3.8) 


(3.9) 
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Furthermore 


(3.10) 


Let 


Then  because  of  F,quation  (3.8) 


dRr  _  rj  -  y 

dr]  ~  Rc 


A  .  d  d 
1  di  J  dll 


V,R 


Rc 


2r:  ~ 


R}RC 


The  surface  element  dA  =  d£zdt]c  transforms  into 
tion  (3.3) 


dC 

M 


dfz  tdR0 

—  1  +  M  a  .  - 
d£  di 


Hence  with  Equation  (3.9) 


(3.11) 

(3.12) 
(3-13) 

d£dr).  One  obtains  from  Equa- 


and  with  Equation  (3.6) 

dJl  -  Rz 

di  -  Rc 

Thus  dA  =  (R0  j Rc)d^dr]  Let 


h{i.nu)  =  h°{£  -  Ut.t),t ) 

Using  Equation  (2.16) 

lU,V,r)  =  *7if  +  jhr,  =  ih°(  +  jh°n  -  5(f\nV) 


(3.14) 


(3.15) 


where  r  is  obtained  from  Equation  (3.1) 

<■  =  <-  ~  =  t- —  me -!)  +  %)  (3.i6) 

and  Rc  is  given  in  Equation  (3.5).  According  to  the  definition  of  h.  Equation  (3.14). 


h(x,yj)  =  h°{x  -  Ut.yJ ) 

then 

hT(x.y.t)  =  h°zo(x  -  Ut.yJ)  =  h°x<,  (x° ,  y°,  t) 
ht(x,yj)  =  -lJhzxo(x  -  Ut.  yJ )  +  hc(  (x  -  Ut.yJ ) 
ht(x.y.t)  =  -Uh:x«(x'  ,y  j)  +  h^{x\y'j) 
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Therefore 

h°t(x'.yc.t )  =  ht(x,y.t)  +  l'hz(x,y,t ) 
One  thus  obtains  the  integral  equation 


u;(a-,7/,0,/)  =  w,:{x.y.0.t)  +  wp(x.y,0.t)  (3-17) 

where 

wp(x,y,0,t)  =  -~(ht(x,y,t)  +  lrhx(x,y,t))  (3.18) 

wv(x,y,0,t)  =  ^  J  I  ~  Mdfi  (319) 

1  =  iht  +  jhn , 


The  vector  p  is  found  in  Equation  (3.7).  When  we  treated  the  problem  in  the  x°,y°,z° 
system,  the  expression  wv  had  been  decomposed  in  order  to  deal  with  the  dominant  sin¬ 
gularity  of  the  integrand  by  analytical  means.  The  same  decomposition  is  carried  out 
here. 

"•  =  =  ■  J  f  mr/V”  <*'20> 

In  the  first  term  a  limiting  process  z  — *  0  is  needed.  According  to  Equation  (3.13) 


P 

RlRc 


=  -v2 


then 


/  /  m?  =  “  / / v*  (i) **•  =  ~  /  <3-21) 

according  to  the  gradient  theorem.  In  the  right  hand  side  of  Equation  (3.21)  one  is  per¬ 
mitted  to  make  the  limiting  process  z  — »  0,  i.e.,  one  replaces  i?0(£,  77 ,  z)  by  R0( £,  77 , 0) .  The 
postponement  of  the  limiting  process  2  — ►  0  has  been  necessary,  because,  in  transforming 
from  £°,  rj°  to  £,  77,  the  contour  of  the  region  will  in  general  not  retain  its  shape;  the  normal 
in  the  £',77°  system  is  different  from  the  normal  in  the  £.77  system  and  so  is  the  length  of 
the  line  element.  The  result,  Equation  (3.21).  is  now  substituted  into  Equation  (3.20). 


wv(x,y,0,t) 


R, 


nd.s  + 


1  f  I  mn.jr)  -  • 7(1,  *,0)^ 
47 T  J  J  R}RC 


(3.22) 
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Note  that  Eq.(3.22)  is  identical  to  equations  (VI. 18)  and  (VI. 19)  of  Reference  1.  The 
transition  from  one  formulation  to  the  other  is  shown  in  Appendix. 

The  effect  of  the  change  of  the  shape  of  an  area  of  integration  connected  with  this 
transformation  is  also  observed  if  one  considers  the  expression  for  the  upwash  in  its  original 
form,  i.i  ..  !>efi>r<  tie  trail-formation-  lo  remove  the  -insularities  are  carried  out.  This 
expres-ion  has  meaning  if  x  and  y  do  not  lie  within  the  region  of  integration.  One  has 
for  the  potential  of  a  source  sheet  in  the  x  .  y  .  £  . //  -system 


os(r  .y  ./)  = 


£  .?/  .T)di  drj 

R 


a 


Hence  in  the  x.  y:  £.  y  systems 


As  before 


oA(.r,  y.t) 


-1  f  r  /?(£.  rj.  T)d£dr) 

4tr  J  I  Rr 


h{i.rt.T )  =  h  (£  -  f'r.y.r) 

T  =  t  ~  al)2  l*W(£  -  X)  +  Rc] 

R<  =  [(£  -  ■r)2  +  -l2(v  -  y V  +  32z2]2 


(3.23) 

(3.24) 

(3.25) 

(3.26) 


Since 


dR(  _  32z 

dz  ~  R, 


one  obtains  in  the  limit  z  —  0 

■.?/.()./)  =  0~:(x .y.Q.t)  =  f  f 


oi(x. 


,2  hU.Tl.T)  hT(Z.T).T) 

Rr  0  Rf 


d^drj  (3d 


Alternatively  one  can  begin  by  forming  the  derivatives  with  respect  to  :  in  the  x  .  y  :  £  .  t; 
svstem-. 


<:>'  !  x  . //  ./) 


l  [  f  h  (£  .y  .-)  M£  .'I  ■-) 


-Jr 


U  ■ 


R1 


d £  d>) 


In  the  subsequent  t ran-forinat ion  to  the  t.-'/  system  one  ha¬ 


h  (£  .  y  ■  ~ )  -  h(i  4  I'm,  .r) 


It. 


For  the  sake  of  clarity  we  denote  the  derivatives  with  respect  to  the  first  and  third  argument 
not  by  he  and  hT  but  by  M1'  and  The  transformation  then  gives  the  alternative 

expression 

hr  =  UhW  +  hW 


therefore 


h{C  rj.r)  r/?(l>(^.r;.r)  +  h(3)(£-  V,  t) 
R2RC  <>  R  RC 


(3.28) 


To  bring  Equations  (3.27)  and  (3.28)  into  agreement  an  integration  by  parts  must  be 
carried  out.  We  are  solely  concerned  with  the  integration  with  respect  to  £,  while  77  is 
kept  constant.  One  notices  that  for  r\  —  constant  is  not  identical  with  M1)  because  £ 
occurs  also  in  r.  We  express  in  the  following  equation  by  means  of  Equations  (3.25) 
and  (3.26). 

=  *i'i  +  hm%  =  hm  _  hmJ 

d£  a0p 2  Rc 

With  this  equation  we  rewrite  the  numerator  of  the  second  term  in  Equation  (3.28) 

uh^  +  *<»>  =  )(i  + 

VhW +hm  =  u£  +  ri^w,s'+ m  ~x)) 


and  with  Equation  (3.6) 


UhM  +h.w  =  u~+h-^ 

d  ^  /if 


Now  Equation  (3.28)  appears  in  the  form 


n  _  1  [  [\h{^n,r)  /i(3)(f,77,7-)l 

02(x,y,O,t)  -  4nJ  j  |  r2Rc  +  R  Rc  +  a  R2  J  dtdr> 


An  integration  by  parts  is  applied  to  the  second  term. 

4(z,y.o.t)  =  ±  [/j^- dv  +  J  JhU.n, 


a  1  h.w 

di  (R0Re)  +  a0R} 


(3.29) 


R  Rc  =  j  _1A/2  [MU-  -r)Rc  + R? 


»  /  1  \  1  1 

d{\RrRj  (l-M2)R2R2 


MRe  +  M^R-)-  +2{£-t) 


17 


Then  one  obtains  for  the  factor  of  h  in  Equation  (3.29) 


1 

R?R? 


+  M-%  1,2  +  2<«  -  x» 


Thus 


,r-A  f2]  .  [»?  +  2-W({  -  *)R'  +  w2u  -  j-)2]  =  a  -  a/2)- 


o- 


f  Mh  ,  [  [ 

'  32h  hT 

’ 

j  rTr/11  +  /  ./ 

hi  +  acR} 

d£di) 

This  agrees  with  Equation  (3.27)  except  for  the  contour  integral.  If  h  =  0  at  the  contour, 
this  integral  vanishes.  The  formula  is  applicable  even  if  h  ^  0  on  the  boundary,  in  which 
case  the  term  accounts  for  the  deformation  of  the  region  which  occurs  as  one  moves  from 
the  £  .  t]  system  to  the  £.77  system. 
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Section  IV 
Lorentz  Coordinates 

As  before  we  denote  x  .  y  .  c  as  a  (.'artesian  system  at  rest  with  respeci  to  the  undis¬ 
turbed  air.  and  with  t  (formerly  /)  the  time.  The  transition  to  wing  coordina ■<  -  he.''  been 
made  by  setting 

x  —  x  -t-  Ala  i  (4.1 ) 

y  =  y 

t  =  1' 

The  Lorentz  transformation"  is  given  by 

x  =  (jc  +  Ma:t:)/3  (4.2) 

y  =  y 

a  t  —  ( a  t  +  M x  )/ 3 

with  the  inverse 

x  =  (x  -  A/a  })/3  (4.3) 

y  =  v 

„c  _  * 

a.i  —  (a/l  —  A  It)/ 3 

It  follows  from  the  first  of  Equations  (4.1)  and  (4.2)  that 

x  =  x/3  (4.4) 

V  =  y 

Sirire  Equations  (4.4 )  do  not  contain  t .  the  j\  y.  c-svstem  is  fixed  with  respect  to  the  x,  y.  z- 
system.  The  time  I  —  t  is  expressed  in  terms  of  t  and  x  by  the  last  of  Equations  (4.3). 
Substituting  here  x  =  ^  one  obtains 

a  t  ~  a  1 3  +  M  ^  (4.3) 
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Equation  (4.4)  and  (4.5)  give  the  relation  between  the  wing  system  j.  y.  z.  t  and  the  Lorentz 
system  x.y.z.t.  One  has  as  inverse  transformation 


x  =  3r 

y  =  y 


[a  t  -  Mi) 
a  '  =  J 

A>  mentioned  above,  rhe  Lorentz  transformation  leaves  the  wave  equation 


d2v  d2x,'  d2v  1  d2v  _ 

dx  2  dy  2  dz  2  a2  dt  2 

unrhanged.  i.e.  one  obtains  in  the  i.y.z.t  system 

d2v  d2v  d2v  1  d2y  _ 

dx 2  +  dy2  +  dz2  a l  dt2 

The  integral  equation  in  the  i.  y,  z,  }  system  therefore  has  the  same  form  as  in  the 
x  .  ;/  .  -  .  t  -system,  namely 


where 


u'(i.y.i) 


o  .i) 


1  dh(x.y.t) 
2°  dt 


R  =  (t  -  x)i  +  (f)  -  y)j 

R  =  R 


~7  =  ih'U.y.r)  +  jhf,(Z,  f/,  r) 

In  this  equation  h(i.y.t)  is  the  unknown  function.  The  transformation  gives 
h(x.y.t)  =  h(x,y.t)  =  h  ^  (ft  a  t  + 

and 

o(x.y.z.t)  =  o(x.y.z.t)  -  0  ^.y.z.  *  (  da  t  -f  j  j 

Hence 

y ,  s,  /)  =  Oj(j.  j/. 0 
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If  o  is  the  velocity  potential,  then  in  most  practical  applications 


<?z(x.y,0.t)  =  u'(x.y.t) 
is  a  given  function,  on  the  wing.  Then 

u'(x.y.i)  —  Qz(x.  y,  0,  /)  =  u-  y  3x.  yja-  ~t  —  Mx)/(al3)  'j 
In  the  indicial  problem  one  has  for  instance 

w(x,y,t)  =  g(x,y)H(t) 
where  the  step  function  H(t)  is  defined  by 


Then 


H(t)  =  I,  t>  0 

H(t)=  0.  1<0 

u'{x,y,t)  =  g((3x,  y)H  (( act  -  Mx)/(a03)) 


Therefore 

u'(x,y,t)  =  g{3x,y).  aj  >  Mx 


■i'(x.y,t )  —  0.  a0t  <  Mx 

The  pressure  difference  between  the  upper  and  lower  side  of  the  wing  is  given  b> 

Ap  =  -p[rhx  +  ht) 


One  has.  from  Equations  (4.4)  and  (4.5) 

/ U \  -  M2  - 

Uhz  +  ht  -  +  —j-  h-t  + 


Therefore 


In  the  wake 


A  P=-%(Uhst  +  hi) 

Uhj  +  h-(  =  0 


The  Lorentz  transformation  maintains  the  general  structure  of  the  integral  equation,  but 
the  expressions  are  somewhat  simpler.  It  requires  a  transformation  of  the  given  up  wash. 


I' or  the  indicial  problem,  the  value  of  /  for  which  ih  jumps  from  zero  to  the  value  g(3x.y) 
now  depends  upon  x.  The  ]>ressures  in  the  x.  y.  t  system  are  expressed  by  the  above  formula. 
After  tin-  pressure-  in  the  x.  'y.  t  system  have  been  found  one  must  transform  back  to  the 
x.  y.l  -ystem.  This  i-  impossible  if  one  tries  to  apply  the  results  in  an  aeroelastic  problem, 
for  there  the  displacements  of  the  wing  are  expressed  in  terms  of  the  time  variable  t  while 
t  he  data  obtained  by  means  of  the  Lorentz  transforms  appear  in  terms  of  t.  This  is  of  minor 
importance,  if  oik  first  evaluates  the  indicial  response.  If  one  uses  the  Lorentz  transform 
then  the  propagation  of  disturbances  occur  with  the  same  speed  in  all  directions.  This 
facilitate-  the  choice  of  the  panels  on  the  wing  and  the  timestep.  especially  at  high  subsonic 
Mach  numbers.  However  the  speed  of  propagation  may  be  of  importance  only  during  the 
initial  phase  of  the  indicial  response.  The  remainder  of  the  report  presents  some  detailed 
discus-ion. 
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Section  V 

Identity  of  Equations  (3.21)  and  Equation  (VI. 18)  of  Reference  1 

The  reader  is  reminded  of  the  following  identity.  Let  £  =  £(l).rj  —  rj{l)  describe  a 
curve  C.  in  the  £r/  plane  and  let  /(£.?;)  be  a  given  one-valued  or  vector-valued  function. 
Then 

^  -  *  ft  Qc  +  d[  drj)dI  -  d!  dl 

In  other  words,  the  operator  t/£  ■  V  represents  differentiation  along  the  curve  C.  Hence 

rlt  -  [B  df 

JA  d*  '  =  JA  dl dl  =  fB  ~ 

If  the  curve  is  closed,  and  /  is  continuous,  then  the  integral  is  zero.  This  holds  for  a  scalar 
as  well  as  a  vector  valued  function: 


(5.1) 


We  address  ourselves  to  the  first  term  on  the  right  of  Equation  (3.20)  after  Equation  (3.21) 
has  been  substituted.  Let 

^  "7  (*.y,0 


and 


4  7T  JR 


I  i  lim  dp 

dz 


(5.2) 


where  is  the  expression  Eq.  (VI.  18)  from  Reference  1.  We  want  to  show  the  identity 
of  /  and  I .  In  the  present  notation,  I  can  be  written  in  the  form 


1  ~  47T1  J  D2n  -  Arr  1  K  X  1 


where 


R?RC 
p(p-dl) 


(5.3) 


(5.4) 


J  f  P\P  ■  d< 

J  R}~R~ 

I  is  rewritten  using  Equation  (3.13) 

1 = -  / ?  K 

To  transform  this  expression,  consider  the  following  integral  which  is  zero  because  of 
Equation  (5.1). 


i 

R 


2?, 


Applying  the  scalar  operator  •  V2  =  +  drj to  the  vector  ^  one  obtains 


''  =  /?  !S{  ( X  )  il  +  Sr,  (  ff  )  H  +  /  i  (§f  ^  +  =  ° 

The  vector  p  is  found  in  Equation(3.7).  One  obtains  this  by  carrying  out  t)  differentiation 
in  the  second  term  and  rewriting  the  first  one. 


MMirW  k 


diT 


+  iM )  d£  +  iM  ~  -  +  j  dr) 


,dRr 


dr) 


=  0 


According  to  Equation  (5.4).  the  first  term  is  —I.  Therefore 

I  -  J  ^-  +  iM  j  [  ~(l°gRo)d(  +  ~{logR0)dr) 

The  second  term  vanishes  because  the  integration  is  carried  out  around  a  closed  curve. 
Therefore 

d* 


=  f 


Ro 


Now 


Thus. 


k  x  d£  =  — fids 

nds 


This  then  shows  the  identity  of  the  expressions  I  (Equation  (5.1))  and  I  (Equation  (5.3)). 
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Section  VI 

Steady  State  Equations 

The  present  formulation  can  be  specialized  to  the  steady  state.  The  equations  so 
obtained  are  different  from  those  obtained  by  a  direct  treatment  of  the  steady  state.  In 
this  section  we  show  by  the  procedure  of  Professor  Williams  that  the  two  formulations  are 
identical. 

In  the  stead.v  state  ir  and  h  are  independent  of  t.  One  obtains  from  Equations  (3.17). 
(3.18)  and  (3.19) 

w(x.  y)  -  wv(x.  ?/)  -t-  wp  (6.1) 


where 


_  -Udh  _  - Mdh 
U  p  2a-  dx  2  dx 


uv  =  li“os/./f  -wi 


F  =  Wit,  =  iU  +  lh 


Substituting  p  from  equation  (3.8)  one  obtains 


1  U  -  t  +  MRC 

h  32  ~~R}RC 

f  __  (*r-  y ) 

2  R2RC 

1  =  ihc  +  jhn 


Written  in  derail  Equation  (C.l)reads 


«•(•*■•  y)  =  ^  Jim/  J  +  f2(Z,r))h  r?]  d£dr)  -  ~hx(x,y) 


Treating  the  steady  case  independently  one  obtains 


i  .  32  ..  f  f  (£  -  r)hi  +  (rt  -  y)hTt 

=  J  - r, i - d(d”  (6  o) 

It  is  our  aim  to  establish  that  the  last  two  equations  are  in  essence  identical.  Because  of 
the  singularity  in  the  denominator,  the  discussion  is  carried  out  for  2  /  0.  Only  at  the  end 
the  limiting  process  2  — >  0  will  be  made.  The  only  procedure  to  transform  one  equation 
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into  tlit*  other  is  to  add  some  expression  to  the  integrand  in  Equation  (6.4)  which  can  be 
transformed,  by  an  integration  by  parts,  into  an  integral  around  the  contour  of  the  region 
of  integration.  (Tins  contour  integral  will  not  always  vanish,  but  at  least  the  area  integrals 
in  Equations  (6.4)  and  (6.5)  can  be  made  to  agree.  The  important  fact  here  is  that  fx  is 
’!'<  factor  wf  he  and  j,  tin-  factor  of  i, ,  I  from  tin  ]>oini  of  view  of  operator  iluori.  one 
makes  use  of  the  relation  that 


j  j  F  ■  (jrad  h  d£d>i  =  contour  term#  -  j  j  h  d.ivF  d£dj/ 

i.e.  that  the  negative  divergence  is  the  adjoint  operator  to  the  gradient.  (This  is  readily 
shown  by  the  integration  by  parts.)  The  double  integral  remains  unchanged,  if  one  adds 
a  vector,  whose  divergence  vanishes,  to  F .  In  the  two-dimensional  case  such  a  vector  can 
be  written  ?  .  where  v  is  the  same  scalar  function.  (In  the  three-dimensional  case, 

it  would  be  curl  A.  where  A  is  some  vector  valued  function.)  Accordingly  we  write,  with 
function  r  (so  far  arbitrary) 

«Y  =  4 \  J  /  [(/i  -  +  (A  +  d^drj  +  f  J -  V'^hjdtdri  (6.6) 

The  second  integral  can  be  transformed  by  integrations  by  part  into  a  contour  integral 
f  / -  Vchn)d(.d})  =  f  h(F  T)){vndri  +  t :(d£)  -  j  j  h{Fn)(Vr,(  -  ip^d^dr) 


The  double  integral  on  the  right  vanishes.  If  /?(£.?/)  vanishes  at  the  contour  of  the  region 
under  consideration,  then  the  contour  integral  vanishes,  too  (otherwise  one  must  carry  this 
expression  along) . 

From  the  postulate  that  the  double  integrals  in  Equations  (6.5)  and  (6.6)  agree,  one 
obtains  as  one  condition 


hit-  v)  +  =  d'1  ( 


v  -  y 

R? 


Hence,  after  substitution  of  f2  from  Equation  (6.3) 


i  f  =  (v  -  y) 


A1  _  1 

/?,’  RcR2 


>i  -  y 


n,R‘  i!R‘  ~  - 


and  after  partial  substitution  of 


(6.7) 


R  -  A  l\M[i  x)  ■+  /?,] 

2U 


(6.8) 


we  obtain 

«•<  =  lM(f  - 112  +  2RM  -  *)  +  MR<1}  <6-9) 

One  remembers  that 

Rc  =  [tf  -  x)2  +  32(t>  -  y)2  +  32z2}  *  (6.10) 

From  this  expression  v  can  be  obtained  by  an  integration;  the  integral  encountered  in  this 
process  exists  in  a  closed  form.  The  simple  result  has  been  given  by  Marc  H.  Williams. 


__  -M{v  -  y ) 

~  RcRc 

This  can  be  verified  by  differentiation.  One  has 


(6.11) 


d{R,Rc)  _  n-2 
di  J 

d(RcRc )  _  1 


R0RC  =  /T2  [Mtf  -  x)Rc  +  Rc2] 
M{i  -  x)2 


MRC  + 


Rr 


+  2(£-x) 


Aircl  =  prc  iMRc 2  +  2R'^  ~x)  +  M(e  -  x>2] 


(6.12) 


_  -1 


[MU  -  x)2  +  2RCU  ~  x)  +  MR2] 


di  32RlR2 
This  vindicates  the  above  choice  of  ip. 

One  still  has  to  examine  the  remaining  terms  in  Equation(6.6)  (and  then  those  in 
Equation(6.4)).  One  has 

to  =  |M(f  _  x)  +  2fl(| 


Then 


^  ~  R0RC  +  {RoRcyRc  +  2Rc ' 


(6.13) 


and  using  Equation  (6.2) 


f\~%  = 

32R.?R2 


__  i  -  x  4-  MRC  M  M(r)  -  y) 


ar-p#(f-x,+“'1 


32RCR2 

fi-4>n=  [(e  -  x)R2  +  MR, ?  +  M32R0R2  -  M02(r,  -  y)2  [Af  (£  -  *)  +  2 Re\] 

Here 

02RO  =  M{i  -  x)  +  Rc 
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and 


32{y  -  y )2  =  Rc  -  u~  x)2  -  32z2 


are  substituted.  Then  the  expression  reduces  to 

r  1  >,9  0wr,  ,  ,  ,  D9l  Mz2  [A/(6  -  x)  +  2RC) 

f]  "  R?32R2  'A/  U  *)  +2MRe(S  x)xRc j  +  £3^2 

The  term  within  the  first  bracket  on  the  right  is  34R2.  Hence 

,  _  32{i-x)  ,  Mz2  [Mtf  ~  x)  +  2RC] 

h  H  '  W  RZR2 

At  least  the  first  term  has  the  expected  form.  So  far  we  have  found 

w{x,y)  =  ^  I  J  ^  -  "—--^3^- — —-didr)  +  j  h(£,  T])(tpr,dri  +  0(d£) 

+ l1™  4nMz  j  j  di  RlR}  L2  +  '  R~  J  d^dr)  2  dx  ^ 


+  l1™  4ttM2  J  J  dj  rIr2  l2  +  j  -  ~2di  (6-14) 

here  0^  and  vv  are  found  in  Equations(6.7),  (6.8),  and  (6.11). 

The  first  term  on  the  right  is  the  one  familiar  from  steady  state  theory.  It  will  be 
shown  that  the  last  two  terms  cancel  in  the  limit  z  — *  0.  The  results  are  summarized  in 
the  form 

4~  j  Jjf \h(_  +  f2hr,)d^dr]  -  ~hz{x,y)  = 

=  |~  /  JA  —  X-ki -  y-htl d£dr)  +  ~  ^  h{i,T)){tpidi  +  ^ydt])  (6.15) 

where  /i ,  /2, 0, 0(,  07,,  i?0?  and  f?c  are  found  in  Equations  (6.2),  (6.3),  (6.9),  (6.7),  (6.11), 
(6.10),  and  (6.9). 

To  show  the  last  two  terms  in  Equation  (6.14)  cancel  in  the  limit  z  —>  0,  one  must 
consider  that  because  of  the  factor  z2,  the  contribution  of  the  second  to  last  term  in 
Equation  (6.14)  vanishes  everywhere  except  for  the  immediate  vicinity  of  the  point  (x,y). 
Thus  can  be  replaced  by|~  and  the  second  to  last  term  in  Equation  (6.14)  simplifies  to 


dh  M 
dx  47t  ' 


(6.16) 


where 


2 

J  =  f  j  r2 r3  \2Rc  +  ~  x-)}d^dy 


(6.17) 
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The  integration  is  extended  over  the  interior  of  a  small  ellipse  given  by 


U  -  -r)2  +  $2{v  -  y)2  -  d2p 


2  -  .k2n2 


pr  —  constant 


(6.18) 


whore  p  ^  0  is  small.  Outside  of  the  ellipse,  the  limit  as  £  — ►  0  is  0.  According  to 

Equations  (3.G) 

R  =  -j2(Rc  +  A/(f-x)) 

Re  =  (U  -  2-)2  +  d2(rj  -  y )2  +  J2r2)i 

(-  X  =  Qz\ 

r)  -  y  =  zy 

2  , 


and  let 


=  2/?(£2  +  r)2  +  1)5 


=  ‘q  (t2  +  f)2  + 1)5 

=  z2  Qd\dr) 

„  2 

Then  the  boundary  of  the  ellipse,  given  by  Equation  (6.18)  becomes  £2  -f  ?)2  =  In  the 
limit  2  — *  0  the  radius  for  this  circle  tends  to  infinity  in  the  \.y- plane.  The  values  of  z  in 
the  integral  cancel  and  one  obtains 


=  //, 


i3{2\/  £2  +  r?2  +  1  +  M£) 


/72 


Now  let 


Then 


and  one  obtains 


where 


(v7e2  +  r)2  +  l  +  +  f)2  +  1)5 

£  =  p  cos  0 
t]  =  p  sin  6 


dt'dy  —  pdpdO 

i  =  d  r  p  i tip) dp 

■I"  (/,2  +  l)-- 

/■  •  *  2 \  />2  -+-  1  +  Af  p  cos  0 

(/»)=/  ' 

7T 


d0 


(\  /»2  +  1  +  MpcosO)2 


(6.19) 


(6.20) 
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The  following  formula  can  be  found  from  the  integral  tables,  for  instance  Reference  3. 


f  A  +  Bcosx  ,  1  r  sin x  [ 

.1  (a  +  b  cost)'*  (T  ~  a2  -  b2  ~  bA)  a  +  bcosx  +  (a'4  "  bB M  a 


dx 

+  b cos x 


It  can.  of  course,  be  verified  by  differentiating  the  right  hand  side.  The  variable  x  cor¬ 
responds  to  0.  For  the  limits  —i\  and  +"  the  first  term  on  the  right  vanishes  and  one 


obtains 


'  T  .4  + .Bcosx  a  A  —  bB  f'n  1 

■  -  -dr=  ai-  b i  .1 ,  i—-' 


r 

J  ,  (a 


(a  4-  b  cos  x)2 


+  b cos  x 


Furthermore  one  has 


r  1 

I  r  o  +  b> 


2  \/a 2  -  b2(tg\) 

- dx  =  ■—= _ =a.rctg - r-  k 

cos  x  Jd/i  -  ft 2  a  +  b 


(again  verified  by  differentiating).  Therefore 


r-n  1 

J-*  a  +  6 CO? 


1  ,  2tt 

—  dx  —  . —  - 

)  COS  X  Ja2  _  f)2 


Hence 


A  +  BcosO  a  A  -  bB 

/  l  r  -  ^dd  =  2n -  3 

J-n  (a  +  b cost))2  f„2_/,2^ 


[a2  -62)s 


Appl  ying  this  formula  to  Equation  (6.20),  one  has 


A  =  2  +  1 


B  =  Mp 

=  \  P2  +  1 
b  =  Mp 


a  A  -bB  =  2  (p2  +  1)  -  M2p2  =  (2  -  M2)p2  +  2 


a2  -  b2  =  p2  +  1  -  M2pl  =  fl2pz  +  1 


f  2  „2  _  /32„2 


Thus 


and.  from  Equation  (6.19) 


I\{p)  =  27r 


(2  -  M2)p2  +  2 
[32  p2  +  1)2 


(6.21) 


7  =  2^r  ((’-‘VV  +  2W 

(/» 


(l>2  +  1)3  (<?V  +  1) 


Lot 


2 

p  —  V 


Th< 
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2  pdp  =  dxi 

*.*  r 

■'0  [  ^,,2 


dll 


(2  -  4-  2 

[  J2  u2  4-  (.L2  4-  1 ) u  4-  l] 2 
which  one  can  find  in  the  tables  of  integrals.  For  A'  =  ax2  4-  2 bx  4-  c 

1  bx  +  c 

A' 2 

1  ax  +  b 
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-- ydx  =  -  ... 

A"  2  or  —  />2 


3  da-  =  -  —  ,, 
5  ac  —  bl 


X2 


In  the  present  case 


a  =  02 

.  $2  +  l  i  M2 

c  =  1 

2  A/4 

ac  —  b  = - -7- 
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With  these  formulae  one  obtains  from  Equation  (6.22) 


_  4t t0  (2  -  2M2  4-  a2/4  -  2  +  2A/2) 
7  ~  M4  3 


I  =  27T 


(6.22) 


To  complete  the  proof  that  the  last  two  terms  in  Equation  (6.14)  cancel,  we  rewrite  Equa¬ 
tion  (6.16)  (the  second  to  last  term  in  Equation  (6.14)). 

dh  M 
()x  4  7T 


and  substitute  our  result  for  I 


which  cancels  with  the  last  term  in  Equation  (6.14). 
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Appendix 

Identity  of  the  Present  Equations  (3.17).  (3.18).  and  (3.22)  with  Reference  1 

In  reference  1  dimensionless  quantities  are  used  while  in  the  present  derivation  the 
quantities  in  their  original  form  are  used.  Equations  (3.17).  (3.18).  and  (3.22)  correspond  to 
Equation  (YJ.22)  in  reference  1.  For  the  present  discussion,  we  characterize  the  quantities 
of  reference  1  hy  a  tilde. 

They  are 


u'Jt.x.y) 


ir('t.x.y)  -  irp(i.x.y)  +  u\.(t .  x.  y) 

-2  7t 


dh  .  -  .  . .  , ,  dh  .  -  .  . . 


W.. 


Here 


(1  -  A/2): 

-  r  +  A/p)/V2)(f,  £.f))  +  (1  -  M2){y  -  y)/t(3>(f .  £.  q)  - 

,  . .  ,?  ‘  "  "  —•  °s ar) 

P  [A/(s  “  j)  +  P)] 

P  =  [(£  -  i)2  +  (r)  -  y)2]  2 


=/r 


(.4.1) 

(.4.2) 


w  — 


do 

dz' 


f  =  t  - 


1 

1  -  A/2 


A/(£  -  r)  +  P 


and  at  the  wing 

o(i.x.y.O)  =  2n'h(i.x,y) 

The  notation  and  h (:{>  refers  to  the  derivatives  of  />(f,£,r})  with  respect  to  its  second 
and  third  argument  respectively.  We  now  turn  to  the  present  formulation.  To  make  the 
analogy  between  the  present  formulation  and  that  of  reference  1  more  evident,  we  write 

Rf  l  0  —  P 

P  =  [u-r)2  +  (l  -A/'Kq-y)2]1 

Then 

R  =  j  \J2  f  v/  (<e  -3 ■)  +  /»] 

The  vector  p  (Equation  (3.18))  then  appears  as 

\2  K  -  T  +  Mp\  +  -  v) 

Furthermore,  we  write  h(r.  £.  r/)  instead  of/?(£.r/.r)  and  denote  by  M2*  and  the  deriva¬ 
tives  of  h  with  respect  to  the  second  and  third  argument.  Then 


V2h  =  ihl2)  +  j7t(;i) 
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Then  Equations  (3.17).  (3.218),  and  (3.19)  assume  the  form 


u'  =  uy  +  tr. 


Up  =  -g-^—  (ht[i.x,y)  +  Uhz(i,  x.  y) 


(.4.3) 


„v  =  1  „..w»)  /  f  <<.  r (.4.4 
4*  ./  J  _  ,q  +  „|12 


/  - 


pIA/(£  -  *)  +  Pj]‘ 

(Af(£  -  x)  +  p] 


a  (l  -  A/2) 

P  =  [(£  -  -r)2  +  (1  -  M2)(r?  -  y) 


2  I  2 


dd> 
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On  the  upper  side  of  the  wing 


0  =  (t,x.y)  =  ^h(t,x,y) 


(see  discussion  preceeding.  Equation  (2.12)). 

The  difference  between  the  two  formulations  is  brought  about  by  the  fact  that  in 
reference  1  the  authors  deal  with  dimensionless  quantities,  and  apply  the  Prandtl-Glauert 
coordinate  transformation.  (In  unsteady  problems  the  Prandtl-Glauert  transformation  is 
less  useful  and  in  a  later  report  Guderley  abandoned  it.)  Moreover  the  definitions  of  h 
differ  by  a  factor  of  4n.  The  quantities  in  Reference  1  are  made  dimensionless  with  a 
characteristic  length,  L,  the  freestream  velocity,  U,  and  the  freestream  velocity  of  sound, 
a. 


One  has 


x  —  xL 

yL 

f]L 

y=  3 

v  =  T 

zL 

~~~  0 

tL 


a, 

0(1,  x,y.  z)  =  U L<t>(i,'x,y.z) 


h(t,x,y)  =  4nULh(i .  x.  y) 
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Therefore 


and 


Since 


Substituting  thes 
(A.l)  and  (A. 2). 


h(t .  x ,  y)  ~  Anxl  Lh 


/t<>  x  y'J\ 
\  L  '  L*  L  ) 


<90 

cte  "  *  ~ 

)/•  =  fTf9?/ 

</£  =  Ldf 

= 

,  Ldii 

rf"=  j 

«Y 

=  Fiiir,. 

1’  =  L  [( f  -  >)2 

+  (;/  - 

v)2]  =  J 

_‘.W,  (■«(*-> 

))  +^j 

1 . 

=  r. 

a 

h(r.  £.  >7)  =  4ttV Lh(r ,  rj) 


/j<2>(r.£.r/)  =  AitVLhW 
hW(T.{.i))  -  4n:iCLh{i) 

(t)  -  t/)M'^(r.  £.  7 7)  =  4 7r r ( f/  -  y)L/<(i)Z, 
(£  —  x  +  A/p)  -  L(l  -  r  +  A/p) 
[A/(£  -  x)  +  p]  =  L  [ A/ ( ^  -  x)  +  p] 


=  ^(2)  =  4jrr/Ai 
ht  =  =  4nUahj 

U  flxvp  =  —2*U(h-t  +  Mh-Z) 

:■  quantities  into  Equations  (A. 3)  and  (A. 4)  one  indeed  obtains  Equations 
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